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We use a relativistic prescription to extract charge and magnetization densities from data for 
the scattering of high-energy electrons by nucleons. A Fourier-Bessel analysis is used to minimize 
the model dependence of the fitted densities. We find that the neutron and proton magnetization 
densities are very similar, but the proton charge density is significantly softer. A useful measurement 
of the neutron charge density is obtained, although the relative uncertainty in the interior will remain 
substantially larger than for the other densities until precise new data at higher Q 2 become available. 
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The electromagnetic structure of nucleons can be in- 
vestigated using the electric and magnetic Sachs form 
factors, Ge and Gm, measured by elastic scattering of 
high-energy electrons. Recent experiments with intense 
high-polarization electron beams have improved the qual- 
ity of the data for nucleon elastic form factors. Perhaps 
the most dramatic observation |l^,[l6| is that the ratio 
between the electric and magnetic form factors for the 
proton decreases sharply for 1 < Q 2 < 6 (GeV/c) 2 . ft 
was suggested that these results demonstrate that the 
proton charge is distributed over a larger volume than 
its magnetization, but radial densities were not obtained. 
At low momentum transfer, the form factors are related 
to Fourier transforms of the charge and magnetization 
densities, but this interpretation is more complicated at 
high momentum transfer. Nevertheless, although models 
of nucleon structure can often calculate the form factors 
directly, it remains desirable to relate the form factors to 
spatial densities because our intuition tends to be based 
more firmly in space than momentum transfer. In this 
paper we use a Fourier-Bessel analysis, together with a 
relativistic relationship between form factors and den- 
sities, to parametrize the nucleon electromagnetic form 
factors in terms of charge and magnetization densities. 

The nucleon electromagnetic vertex function takes the 
form 
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where F\ and F 2 are known as Dirac and Pauli form 
factors, K is the anomalous part of the magnetic mo- 
ment, and 7 M and a^ v are the usual Dirac matrices {e.g., 
0). This current operator appears simplest in the Breit 
frame where the nucleon approaches with initial momen- 
tum — <7b/2, receives three-momentum transfer qg, and 
leaves with final momentum (fs/2 such that the energy 
transfer vanishes. In the Breit frame for a particular 
value of Q 2 , the current separates into electric and mag- 
netic contributions [HI 
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where Xs is a two-component Pauli spinor and where the 
Sachs form factors are given by 



G E = Fi- tkF 2 
G m =Fi + kF 2 
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with r = {Q/2m) 2 . Naively it would appear that 
one could obtain charge and magnetization densities as 
Fourier transforms of the electric and magnetic form fac- 
tors, but each momentum transfer Q specifies a different 
Breit frame. Early experiments with modest Q 2 
gested that 
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where \i is the magnetic moment in nuclear magnetons 
and G D {Q 2 ) = (1 + Q 2 /A 2 )- 2 with A 2 = 0.71 (GeV/c) 2 
is known as the dipole form factor. Similarly, data 
for Geu at low Q 2 can be described by the Galster 
parametrization |3j 



G En {Q 2 ) ^ -p n G D {Q 2 ) 



1 + Bt 



(4) 



where A and B are constants. 

Let p c h{r) and p m {r) represent spherical intrinsic 
charge and magnetization densities normalized according 
to 



dr r 2 p ch (r) 



dr r 2 p m {r) = 1 



where Z = 0, 1 is the nucleon charge. An intrinsic form 
factor can then be defined by the Fourier-Bessel trans- 
form 



P{k) 



dr r 2 j (kr)p{r) 
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where k is the intrinsic spatial frequency (or wave num- 
ber). If one knew how to obtain p{k) from data for the 
appropriate Sachs form factor, the intrinsic density could 
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be obtained simply by inverting the Fourier transform, 
such that 



p(r) 



dk k jo{kr)p(k) 
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The naive or nonrelativistic inversion from form factor to 
density replaces Q by k and p(k) by G(Q 2 ), but leads to 
unsatisfactory densities for dipole or Galster form factors 
- the direct Fourier transform of a dipole form factor 
is an exponential density with a cusp at the origin. The 
corresponding density for the Galster form factor is more 
complicated but also features an unphysical cusp at the 
origin. 

Although there is no rigorous model-independent re- 
lationship between the Sachs form factor and the cor- 
responding static density in the rest frame, a variety of 
models offer prescriptions of the form 
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where A is a model-dependent constant and where the 
intrinsic spatial frequency k is related to the invariant 
momentum transfer by 



k 2 = 



Q 2 
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The most important relativistic effect is the Lorentz con- 
traction of spatial distributions in the Breit frame and the 
corresponding increase of spatial frequency represented 
by the factor of (1 + r) in Eq. (||). A measurement 
with Brcit-frame momentum transfer qs = Q probes 
a reduced spatial frequency k in the rest frame. The 
Sachs form factor for a large invariant momentum trans- 
fer Q 2 is determined by a much smaller spatial frequency 
k 2 = Q 2 /(1 + t) and thus declines much less rapidly with 
respect to Q 2 than the Fourier transform of the density 
declines with respect to k 2 . Licht and Pagnamenta Q] 
demonstrated that by accounting for Lorentz contraction 
a good fit to data for Ge p can obtained using a Gaussian 
density typical of quark models. Thus, this procedure 
provides physically reasonable densities free of cusps at 
the origin. However, according to Eq. (||), the maximum 
spatial frequency that is accessible with spacelike Q 2 is 
k m = 2m. The limitation k < k m can be interpreted as 
a consequence of relativistic position fluctuations, known 
as zitterbewegung, which obscure details finer that the 
nucleon Compton wavelength. 

The first attempt to relate elastic form factors to 
ground-state densities was made by Licht and Pagna- 
menta [pi using a quark cluster model in the impulse 
approximation and a kinematic boost. They proposed a 
relativistic inversion method using Eq. (j^) with Xe = 
\m = 1- However, these choices for A do not ensure that 
the Sachs form factors scale with Q~ 4 , as expected from 



perturbative QCD, unless restrictions are placed upon 
p(k m ). Using a more symmetric treatment of the cluster 
model, Mitra and Kumari || found that A^ = Am = 2 
automatically satisfy the perturbative QCD scaling rela- 
tions at very large Q 2 without constraining p(k m ). More 
recently, Ji obtained similar relationships with Xe = 
and Xm = 1 using a relativistic Skyrmion model based 
upon a Lorentz invariant Lagrangian density for which 
the classical soliton solution can be evaluated in any 
frame. Quantum fluctuations were then evaluated after 
the boost. Furthermore, the recent soliton calculation by 
Holzwarth (?]] uses this prescription and accurately pre- 
dicted the decline in Ge p /Gm p - In this paper we employ 
the values Xe = and Am = 1 suggested by the soliton 
model and in a subsequent paper || will investigate in 
detail the discrete ambiguity due to the choice of A. 

To extract radial densities from the nucleon form fac- 
tor data we employ techniques originally developed for 
fitting radial distributions to data for scattering of elec- 
trons or protons from nuclei Js| — jXlj] - Simple models with 
a small number of parameters do not offer sufficient flex- 
ibility to provide a realistic estimate of the uncertainty 
in a radial density. Rather, we employ linear expansions 
in complete sets of basis functions that are capable of de- 
scribing any plausible radial distribution without strong 
a priori constraints upon its shape. Such expansions per- 
mit one to estimate the uncertainties in the fitted density 
due to both the statistical quality of the data and the in- 
evitable limitation of experimental data to a frequency 
range, k < k max . The uncertainty due to limitation of 
k is known as incompleteness error. More detailed dis- 
cussion of the method may be found in Refs. g-[yj|, but 
the basic idea is to supplement the experimental data 
by pseudodata of the form p(ki) = ± Sp{ki) whose 
uncertainties are based upon a reasonable model of the 
asymptotic behavior of the form factor for fej > k max 
where k max is the spatial frequency corresponding to the 
maximum measured Q 2 . On quite general grounds one 
expects the asymptotic form factor for a confined system 
to decrease more rapidly than fc~ 4 []l0| . Therefore, we 
assume that 



5p(k) = p(k max ) 



(9) 



Notice that we must apply this procedure to p(k) instead 
of the Sachs form factors because, according to Eq. (^), 
the accessible spatial frequencies are limited to k < 2m 
for spacelike Q 2 while convergence of Eq. (^) requires a 
constraint upon p(k) for inaccessible range k > 2m. 
The Fourier-Bessel expansion (FBE) takes the form 



P( r ) = ^2a n j (k n r)Q(B, - r) 



(10) 



where is the unit step function, R is the expansion 
radius, k n = nir/R are the roots of the Bessel function, 
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and a n are the coefficients to be fitted to data. One ad- 
vantage of the FBE is that the contribution of each term 
to the form factor is concentrated around its k n so that a 
coefficient a n is largely determined by data with k ~ k n . 
The larger the expansion radius R, the smaller the spac- 
ing between successive k n and the greater the sensitivity 
one has to variations in the form factor. One should 
choose R to be several times the root-mean-square ra- 
dius but not so large that an excessive number of terms 
is needed to span the experimental range of momentum 
transfer. Terms with k n > k max provide an estimate of 
the incompleteness error. We chose R — 4.0 fm, but 
the results are insensitive to its exact value. Small but 
undesirable oscillations in fitted densities at large radius 
were suppressed using a tail bias based upon the method 
discussed in Ref. |L2]. We employed a tail function of 
the form t(r) oc e~ , based upon the successful dipolc 
parametrization for low Q 2 , and included in the x 2 fit a 
penalty for strong deviations from the tail function for 
r > 2.0 fm. The constraint on the neutron charge was 
also enforced using a penalty function. The tail bias im- 
proves the convergence of moments of the density but has 
practically no effect upon a fitted density in the region 
where it is large. The error band for a fitted density is 
computed from the covariance matrix for the \ 2 fit an( i 
includes the incompleteness error, but does not include 
discrete ambiguities due to the choices for A in Eq. (]?]). 

We selected the best available data in each range of 
Q 2 , with an emphasis upon recent data using recoil or 
target polarization wherever available. Gm p data were 
taken from the compilation of Holder [ffsl for Q 2 < 0.15 
(GeV/c) 2 and for larger Q 2 from the analysis of Brash 
et al. 1 14 1 using the recent recoil polarization data for 
Gev/Gmt, from Refs. ]l5||L6[ |. Cross section data from 
Refs. were used for Q 2 < 1 (GeV/c) 2 but cross 

section data for Ge p were excluded for larger Q 2 . Simi- 
larly, the data for Geu were limited to recent polarization 
data ||l9|-p4|, with corrections for nuclear structure and 
final-state interactions whenever available, plus the anal- 
ysis of and T20 by Schiavilla and Sick pq| . We also 
included the neutron charge radius obtained by Kopecky 
et al. |5q| from transmission of thermal neutron through 
liquid 208 Pb and 209 Bi. Finally, for Gmu we selected 
polarization data from pTfl and cross section data from 



Fits to the form factor data are shown in Fig. [l] as 
bands that represent the uncertainties in the fitted form 
factors. The widths of these bands are comparable to 
the experimental precision where data are available, but 
expand for larger Q 2 where the uncertainties are based 
upon Eq. (||). Notice, however, that in this figure the 
form factors were divided by Gd, which decreases rapidly 
with Q 2 . Therefore, the impact of uncertainties at large 
Q 2 upon the fitted densities remains modest because, 
with the exception of Ge„, the form factors themselves 
and their absolute uncertainties become rather small at 



their largest Q 2 . Although the low-Q 2 data for Gmu 
have improved in recent years, significant systematic dis- 
crepancies remain. Recent data from Refs. p?| , |3l| , p3| , ^| 
with small statistical uncertainties suggest a small dip 
near 0.2 and a peak near 1 (GeV/c) 2 . For Geh we plot 
the Galster model Q for comparison. The simple two- 
parameter fit Galster et al. made to the rather poor data 
available at that time did not permit a realistic estimate 
of the uncertainty in the form factor or fitted density and 
the apparent agreement with more modern data must be 
judged as remarkable but fortuitous. 

Proton charge and magnetization densities are com- 
pared in Fig. Both densities are measured very pre- 
cisely. The new recoil-polarization data for Ge p decrease 
more rapidly than either the dipole form factor or the 
magnetic form factor for Q 2 > 1 (GeV/c) 2 ; consequently, 
the charge density is significantly softer than the mag- 
netization density of the proton. Neutron densities are 
shown in Fig. We find that the magnetization den- 
sity for the neutron is very similar to that for the proton, 
although the interior precision is not as good because 
the range of Q 2 is smaller and the experimental uncer- 
tainties larger. Limitations in the range and quality of 
the Geu data presently available result in a substantially 
wider error band for the neutron charge density. Data at 
higher Q 2 are needed to improve the interior precision, 
but a useful measurement of the interior charge density 
is obtained nonetheless. The positive interior density is 
balanced by a negative surface lobe. Note that polariza- 
tion measurements are sensitive to the sign of the density, 
but that cross section measurements are not. 

In summary, we have developed a parametrization of 
nucleon electromagnetic form factors based upon radial 
densities expanded in complete bases. Although we can- 
not claim that the fitted densities are unique, this method 
does provide an intuitively appealing parametrization of 
the data. It is necessary to account for Lorentz con- 
traction at high momentum transfer in order to obtain 
reasonable interior densities without a cusp at the ori- 
gin. Using recent recoil polarization data for the proton, 
we find that the charge density is significantly broader 
than the magnetization density. We also find that the 
neutron and proton magnetization densities are similar. 
Although the available data are not yet as precise, the 
neutron charge density is also obtained with useful preci- 
sion. A more detailed discussion of the model dependence 
of fitted densities will appear elsewhere. 
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